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ABSTRACT. A discussion of the FSRQ - BL Lac unification. All distinctive features marking 
these two Blazar subclasses find an unifying explanation if the sources are powered by central 
engines constituted by similar Kerr holes, but fueled at high and low accretion rates, respec¬ 
tively. The connection need not be a genetic one, but evidence toward some FSRQs switching 
into BL Lacs at lower 2 will be provided by moderately negative BL Lac evolution. Then an 
extrapolation will be warranted toward ultra-high energy particle accelerators operating at very 
low accretion rates. 


1. Introduction 

The BL Lac Objects may be described as minimal Blazars. This is because - unlike the 
FSRQs which constitute the other main Blazar subclass - they feature at the same time 
weak or no emission lines and blue-UV bumps, weak or moderate intrinsic power, and 
little or no cosmological evolution (for general descriptions and detailed references, see 
Urry & Padovani 1995, Padovani & Urry 2001). Being so rare, elusive, with redshifts 
not easily pinpointed, and hence so vulnerable to selection effects, these sources make 
difficult observational targets. 

Yet the BL Lacs constitute the focus of keen and mounting interest, being extreme 
Blazars that shine by nearly pure non-thermal radiation with little or no photon repro¬ 
cessing. Their spectra are exceedingly wide, in that they cover 16 frequency decades or 
more, from the radio band (as also the FSRQs do) up to the high energy 7 -rays beyond 
10 GeV. In fact, a number of such objects have been observed into the TeV range (see 
Costamante et al. 2002). 

Concerning the spectral energy distribution (SED), all Blazars follow a roughly sim¬ 
ilar pattern with two bumps mainly or mostly contributed by Synchrotron and Inverse 
Compton emissions. But (as pointed out by Fossati et al. 1998) in the FSRQs the fre¬ 
quency Vp of the lower peak is placed in the bands from the far IR to the optical, and 
the upper peak is around 10^ MeV. In the BL Lacs, instead, Vp ranges from optical to 
X-rays, and the upper peak (possibly modified by the Klein-Nishina limit) falls around 
or above 10 GeV. 

In all Blazars such non-thermal emissions are sharply beamed; they arise from rel¬ 
ativistic jets of particles with bulk Lorentz factor F around 10 (see Sikora 2001), not 
substantially different for BL Lacs and FSRQs. These jets are seen nearly pole on within 
angles of order l/F, with apparent powers enhanced by a “blazing” factor F^ relative to 
the emitted ones (see Sikora et al. 1997, Ghisellini 1999). What differs between the two 
Blazar subclasses is the total output produced; the top outputs are apparently limited 
to L ~ 10^® erg/s in the BL Lacs, but may be up to 10^ times stronger in some FSRQs 
(see Maraschi 2001). 

How and why all these BL Lac peculiarities arise together? Largely on model-indepen- 
dent grounds, we trace back their common origin to different fueling levels of basically 
similar power sources. 

2. The engine 

The paradigm for the energetic cores of all AGNs is of course provided by accreting 
black holes with masses M ~ 10®^^ Mq. The accreting gas gathers in a disk where its 
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angular momentum is transferred outwards; in the process thermal, roughly isotropic 
emissions are generated with power Lth ~ rjLErn in terms of the conversion efficiency 
r] ~ 10“^, of the Eddington luminosity Le and of the corresponding accretion rate rfi. 
Does the latter differ markedfy in BL Lacs and in FSRQs? 

An affirmative answer is provided evidence concerning the thermaf emissions. The 
FSRQs, fike the radioquiet QSs, emit a powerful “big blue bump” conceivably from 
the inner disk, in addition to strong lines from continuum reprocessed by distant gas 
“clouds”; the luminosities are high (Tavecchio et al. 2000) and close to the Eddigton 
values for BH masses around 10® Mq, so to ~ 1 — 10 is to hold. 

In BL Lacs, instead, the absence or weakness of both the emission lines {EW ^ 5 
A) and the BBB call for conditions of low particle densities both around and in the 
accretion disks. Since the latter density scales as n oc to^^/®® and the former as n oc to 
( see Frank, King & Raine 2002), these conditions indicate to ^ 1. A similar conclusion 
holds if the BBB is due to emission from a hot corona (Sun & Malkan 1989). 

On the other hand, evidence from the jets and their emission indicates the mass is not 
the only parameter relevant to the central BH in Blazars. Such an evidence is provided, 
e.g., by the strong radio emission powered by jets which has been recently shown not 
to be simply correlated with large host and BH masses (Ho 2002, McLure & Dunlop 
2002). So a further variable must enter; we entertain the view involving the other basic 
parameter of the BHs, i.e., the angular momentum J (see Blandford 1990, 1993). 

By its vectorial nature this is very likely to provide the sharp and steady directionality 
of the jets. In addition, rotating Kerr holes make closer stable orbits (approaching Vg « 
GM/c^ w l.SlO^'^Mg cm) available to the gas particles before they plunge into the 
horizon; so more gravitational energy can be extracted, and the the maximal efficiency 
is raised from 77 ~ 6 10“® to nearly 0.4. Finally, the rotational energy of a Kerr hole 
(proportional to J® in a classical rendition) may contribute directly the jet-like output 
of BL Lacs in particular. 

Part of the rotational energy of the hole can be directly extracted via hydromagnetic 
Poynting-like flux related to B and E fields coherent on large scales, in approximately 
force-free conditions E • B w 0; this is of course the attractive but debated scheme 
proposed by Blandford & Znajek (1977), with the power produced Lk ultimately de¬ 
pending on the square of the B field rooted in the disk but threading the horizon close to 
Tg. A similar hydromagnetic mode of energy extraction has been extended by Blandford 
& Payne (1982) to the contribution Ld from the accretion disk itself. 

Thus two contributions to the non-thermal power Lk + Ld = Lnth may be envisaged 
from the central region surrounding the Kerr hole within radii re ^ Cg, in the presence 
of the poloidal field B. Both components scale approximately as 

Lk cx B'^rl , ( 1 ) 

following the basic dependences provided by electrodynamics. But the pressure P = 
B'^ I8tt holding the magnetic field, and the specific value of the effective radius Ve are 
provided by the disk structure, and this depends on the accretion rate to. 

3. Fueling rates 

So both outputs Lnth and Lth are ultimately governed by the fueling regime of similar 
basic engines. But the latter output scales down as Lth oc to when to is decreased, while 
the same may not be the case with the former. 

Dependences and numerical coefficients in Lnth have been focused and discussed 
by Ghosh & Abramowicz (1997); Moderski, Sikora & Lasota (1998); Livio, Ogilvie & 
Pringle (1999). In the framework of standard disk models the BZ77 scheme is found to 
yield the maximal power 


Lk^2 10^^ Mg { J/Jm? erg 


(2) 
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where Jm ~ GM^/c is the maximum value for the angular momentum consistent with 
the Kerr metric. This value of Lk is attained for accretion rates m ^ 10“^ when the 
inner disk is dominated by the radiation pressure P oc M~^ independent of rh. 

For lower m instead, the gas pressure P oc dominates also the inner 

disk; then lower powers Lk « 10'^'* {J/ JmY erg/s obtain. Yet the the ratio 

Lx/Lth exceeds unity and scales inversely with to, as to“^/®. The scaling is stronger if 
indeed at to w 10“^ a silent ADAF regime sets in (see Narayan, Mahadevan & Quataert 
1998; Frank et al. 2002); then the ions cannot share much energy with electrons before 
plunging into the horizon, and most radiations are suppressed sharply. 

We stress two main points of interest in these results. First, the emission becomes 
increasingly non-thermal anyway as to decreases below some 10“^. Second, for values 
TO > 10“^ the power saturates to Lk ~ 6 10^® erg/s for M ^ 3 10® Mq. To this value 
one should add the the nonthermal disk contribution Ld, which in conditions of maximal 
J and closest orbits with small Tg « rg may be down to the same order as, or just a few 
times larger than Lk (Meier 2001). Then the Kerr hole and inner disk may be considered 
as a dynamically and magnetically coupled system, using up also energy stockpiled in 
previous accretion episodes of mass and associated angular momentum (Bardeen 1970). 

Thus a total Lnth ~ 10^® erg/s appears to provide an upper bound to outputs mainly 
constituted by non-thermal radiation, the defining feature of a BL Lac; the condition 
Lnth/Lth ~ 1 remains relevant in spite of the blazing effects of the jets that tend to 
swamp a comparable thermal emission. This leads us to propose that top BL Lac outputs 
remaining below this limit provide evidence that rotational energy is in fact extracted 
from the hole via the BZ77 mechanism. The existing data (Maraschi 2001, see also fig. 
1 ) support to now this prediction. 

The FSRQs, on the other hand, do have larger outputs that in some cases exceed 10^^ 
erg/s if protons contribute more then the electrons to the beam kinetic power (Celotti, 
Ghisellini & Padovani 1997; Tavecchio et al. 2000). Within the above framework, such 
outputs require a dominant contribution Ld 3> Lk from a wider disk region dominated 
by radiation pressure; this in turn requires conditions where to ~ 1 , consistent with the 
independent requirement posed by the high levels of Lth discussed in §2. 

The maximal power expected may now approach 10^® erg/s, exceeding the maximal 
Lk by some 10® (Livio et al. 1999). This obtains if the magnetic field in the disk 
is and remains a few times larger than the value threading the hole horizon out to 
distances r^ w 5rg. Here both Ld and the thermal emission Lth are directly fed by 
current accretion; both scale as to, and Lnth/Lth ~ 1 is expected to hold. Now the hole 
contribution is subdominant, but is likely to provide a “high-velocity spine” crucial for 
the outward propagation of the jet (see Livio 1999, Chiaberge et al. 2000). 

The detailed share between Ld and Lth depends on the balance between two scales 
in the power spectrum of the magnetic field inhomogeneities: the large-scale, coherent 
vs. the small-scale, turbulent component. A dominant share of the former, such as to 
yield a large Ld and correspondingly an additional, large transfer of angular momentum 
outwards, may require modifications of the standard disk models (see Salvati 1997). 

Alternatively, to account for the huge FSRQ outputs one would need very strong B 
threading the hole, up to the values B"^/Stt ~ pc® in the plunging orbit region advocated 
by Meier (1999). Such enhanced fields have been argued and discussed variously; they 
look unlikely in a thin disk, whereas in thick disks their status is still uncertain. We 
note they would require high levels of to anyway, the main issue here. 

In sum, independently of the open if interesting questions concerning the detailed 
modes of power production, these considerations strongly suggest a trend toward weaker 
powers L but higher ratios Lnth/Lth as to is decreased from values 1-10 that mark the 
FSRQs to 10“® that mark the BL Lacs. We (CD02) referred to this trend as the Blazar 
Main Sequence (BMS), with the implication that it holds independently of any genetic 
link between BL Lacs and FSRQs. 
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Fig. 1. The L vs. Vp diagram. Data from Costamante et al. (2002). Even the most powerful BL 
Lacs measured to now fall below L ~ 10“^® erg/s. The solid line represents the trend L oc Vp^^^ 
provided by eq. (3). The dashed line represents the intrinsic scatter predicted in §4a. Crosses 
indicate the positions of Mkn 421 and Mkn 501 at the top of their flares. The recent data by 
Giommi et al. (2002a) add substantial scatter particularly in the BL Lac section of the diagram. 


4. The Blazar Main Sequence 

The BMS may be usefully contrasted with the stellar MS. Stars were placed, even 
before their energy source was properly understood, on the HR diagram relating the 
simplest observables: luminosity and color, or black body peak izi, (x Tg. Now the MS 
comprises the sources whose luminosity is powered by central nuclear burning of H. In 
stars, the spectra of the outgoing radiation are eventually shaped in the upper envelope 
under conditions still close to optically thick. All processes occur in conditions close to 
hydrostatic and thermal equilibrium; so the spectra are closely black body, looking like 
spikes at Vb on a broad frequency span, while L increases sharply (and the lifetimes 
decrease) with increasing i^p. The main underlying parameter is mass. 

Blazars, instead, are ruled by strong gravity; in addition, the jet-like, nonthermal 
component of their radiation is emitted by highly relativistic particles under optically 
thin conditions. The lower peak at Vp in their SEDs is broad, yet sufficiently well de¬ 
fined to be represented on the plane labeled by the simplest observables: Vp itself, and 
the integrated steady luminosity L associated with the peak. Thus a definite “spectral 
sequence” has been obtained by Fossati et al. (1998), with the powerful FSRQs at low 
Vp and the weak BL Lacs at high Vp] the latter data have been recently extended by 
Costamante et al. (2002). We hold to to be the main parameter marking the BL Lacs 
from the FSRQs (see also Bottcher & Dermer 2002), and underlying the whole BMS. 

We stress analogies and differences between the BMS and the stellar MS on represent¬ 
ing the Blazar data first in the simplified form of fig. 1. The key difference is constituted 
by the decreasing trend of L vs. Vp in the BMS as compared with the increasing trend 
in the stellar MS; this BMS feature is resistant to selection effects and new detections 
presented and discussed by Giommi et al. (2002a). We trace back such a trend to con¬ 
ditions close to steady state (steady accretion and steady jet outflow) underlying the 
BMS, that replace the equilibrium conditions governing the MS. 

In fact, as a condition for steady energy distribution N{j) of the relativistic electrons 
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that emit S and IC radiation, one may equate the acceleration and the radiative cooling 
times, taking up from Ghisellini et al. (1998) and Bottcher & Dermer (2002). So from 
7 /i? oc 1/jL one obtains the decreasing course L oc at given accelerating electric 
field B. 

But we find a flatter dependence when the primary B fields are stronger, due to 
lower densities n oc rh in the source. For example, consider the parallel fields E • 
B 0 which arise at the flow boundaries where breakdown occurs for the force-free 
approximation governing the bulk of the BZ77 magnetospheres; this is bound to occur 
in a space-inhomogeneous and also in a time-depending fashion on short scales. Then 
electrodynamic screening produces effective energy gains B d (x c/ujp oc that scale 

up when the electron density in the plasma frequency wp is lower; so the steady state 
constraint oc I/ 7 T will allow larger 7 for a given L. 

To see how the course of L vs. is modified, one may iteratively relate the minimum 
density to emissivity to obtain n <x L\ matters are simplified on keeping the source 
volumes roughly constant. Now the course flattens to 

L oc , (3) 

not a bad rendition even for the overall trend in the data, see fig. 1. To a better ap¬ 
proximation, one may take into account (see CD02) the full dependencies of tOp, n,and 
L/n (X on the minimal or maximal electron energies 7 ^, 7m and on the slope of 
to obtain a sharper flattening L(xvp in the BL Lac range. Such flatter courses may 
be also viewed at as “weaker” correlations. 

By the same token, we expect the maximal electron energies to scale up following 
7m ^ eB d/nie (? ~ 10® i ?4 dio (Mg/riy)^-^® oc -7 here the maximum field is 

expressed as E ^ 10 ^i ?4 G and (following Blandford & Payne 1982) is taken to scale 
as E oc into the emission region at r ~ 10 ^®^^ cm, while the effective screening 

distance is d = 10^® dio cm. So we obtain values 7 m ~ 10® in FSRQs, and values up 
to 10^ times larger in BL Lacs; these are high enough to produce by IG radiation the 
observed TeV photons. 

We add that N{'j) is expected to depart from a pure power-law as the acceleration 
takes place in fields endowed with some degree of coherence. In fact, we expect the 
electrons to be accelerated on crossing a sequence of many sheets or filaments where the 
force-free condition gradually breaks down on average. Gorrespondingly, the parallel E 
field will progressively grow to a point while the width of the filaments or sheets widens, 
with a corresponding growth of the energy gain per acceleration step. 

On the other hand, conditions where the energy gain exceeds a factor 2 (a limiting 
value in the acceleration by relativistic shocks, see Achterberg et al. 2001) to approach 
10 are discussed by Massaro (2002). Then the energy distribution N{j) has a eurved 
shape, actually a log-parabolic one; the emitted S or IG spectra are correspondingly 
curved over a wide frequency band. In fact, the data concerning a number of BL Lacs 
apparently show such a curvature in the O - X-ray range (Massaro 2002, Giommi et 
al. 2002 b); so they call for energy gains of order 10 , more appropriate for accelerating 
fields with some degree of coherence. 

Having outlined the simplified picture, next we point out a number of additions 
yielding variance, a) The BMS has an intrinsic width related to the additional parameter 
LxlLd that implies within the BL Lac subclass some asymmetric scatter toward high 
L. This is because only Ld is directly related to m and to n and (inversely) to Vp by 
the arguments in §3, while Lk tends to be independent when relevant at low to; then 
higher non-thermal luminosities Lnth = Ld + L^ obtain in BL Lacs at a given Vp. b) We 
recall that we expect 7m oc to hold in the BL Lac range; here we add that B is 

likely to be distributed among the objects with a bias toward low values. This will yield 
lower values of Vp and L, contributing more asymmetric scatter, c) Additional scatter 
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Fig. 2. The BL Lac counts evaluated from the expression given by CD02. The short-dashed line 
represents the result from pure LE with time scale tl = 7 Gyr. The solid line represents the 
result when negative DE with scale td = —5 Gyr is included. We show also the ESRQ counts 
from the same expression, with tl ~ 2 Gyr and td = + 5 Gyr. The dotted line represents 
the “Euclidean” slope. Hubble constant Ho = 65 km/s Mpc. Data for BL Lacs from Giommi, 
Menna, & Padovani (1999). 


toward low L is generally contributed by misaligned sources with a given value of T (see 
Georganopoulos 2000), and by the distribution of F itself. 

These different scatter components will superpose in the BL Lac range to a flattening 
“bare” correlation, to blur it mainly toward low values of L and Vp in a manner not 
unlike that observed and discussed by Giommi et al. (2002a). 

Finally, episodes of violent variability are observed to occur in the form of flares from 
hours to weeks, particularly in BL Lacs at high photon energies (see Costamante et al. 
2002). Here the steady state constraint clearly does not apply; rather, the basic transient 
behavior of the S and IC radiations L oc 7 ^ oc t'p provides a first approximation to the 
departing branches (see fig. 1) traced by flaring objects on the L — Vp plane. 

5. Cosmological evolutions 

Along the BMS, we expect longer object lifetimes and slower cosmological evolution to 
occur. Basically this is because conditions of weak activity feeding on low levels of m 
can live longer than those of hyperactivity requiring high to. 

This may be quantified on considering that roughly similar BH masses of order 
M ^ 3 10® Mq found (see Treves et al. 2002) in BL Lacs and in FSRQs imply a similar 
limit to cumulative masses accreted during the activity of either type. As to the latter, 
the mass constraint bounds the cumulative activity time to a few Gyrs with a duty cycle 
of order 10“^ (Gavaliere & Padovani 1989). 

As to the former, instead, the constraint allows a lifetime up to 10 Gyr at levels 
lO"*® erg/s or weaker; note that a duration E/Lk ~ 5 Gyr can be sustained just by 
the rotational energy extractable from the coupled system Kerr hole - inner disk, in 
the absence of angular momentum replenishment (Gavaliere & Malquori 1999). Since 
a similar time behavior L{t) is shared by these objects at various power levels, such 
values apply to the specific scale tl ~ 10 Gyr for the “luminosity evolution” of their 
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population. 

On these simple grounds we expect the space distribution of the BL Lacs to be uni¬ 
form to a first approximation, out to some z « 1; so we expect from the volume test val¬ 
ues (K/Oa) « 0.5, and for the number counts a nearly Euclidean shape N{> S) oc S~^'^ 
at intermediate fluxes. This, in fact, what the data indicate to a first approximation, 
see Giommi, Menna & Padovani (1999). 

A quantitative prediction is obtained by CD02 on using the expression for the counts 
at intermediate fluxes in terms of evolutionary times, in turn evaluated in more detail. 
The results for BL Lacs counted in the radio band are plotted in fig. 2; the short-dashed 
line is obtained on using only the basic time scale tl ~7 Gyr that marks their slow LE. 

In contrast (see the long-dashed line in fig. 2), quite steeper counts obtain for the 
FSRQs using the time scale w 2 Gyr appropriate for their strong evolution, similar 
to that of radioquiet QSs (Goldschmidt et al. 1999; Giommi, Menna & Padovani 1999). 
The latter behavior is widely traced back to exhaustion on the scale of a few Gyrs of 
the gas stockpiled in the host galaxy and usable for accretion; this is caused by previous 
accretion episodes, in addition to ongoing star formation (Gattaneo, Haehnelt & Rees 
1999; Haehnelt & Kauffmann 2000; Gavaliere & Vittorini 2002). Exhaustion is bound 
to occur at z < 2.5, when the formation of cosmic structures by hierarchical clustering 
evolves beyond the galactic scales; then violent merging events assemble the galaxies 
into richer and richer groups, but rarely reshuffle the galactic masses and no longer 
import much fresh gas supplies into the host galaxy. 

In fig. 2, another time scale scale appears, namely, td « ± 5 Gyr; this is negative for 
the BL Lacs (solid line), meaning negative “density evolution”, that is, birth of these 
objects at moderate or low z. In these conditions, less and less sources are found on 
looking out to increasing z; this implies 14/14 < 0.5 and flattening of the counts below 
the Euclidean slope at relatively high fluxes, where the low redshifts of the BL Lacs still 
prevent the cosmological convergence of the volumes from operating. 

Negative evolution is not yet apparent in all existing the BL Lac surveys (see Cac- 
cianiga et al. 2001, Perri et al. 2002), but it will eventually provide an interesting telltale 
about the origin of these peculiar AGNs. To illustrate the issue, let us consider a genetic 
link between the Blazar subclasses such that BL Lac birthrate w — FSRQ deathrate. 

First, consider that such a link is a conceivable one, since the basic engines are similar 
while the fueling rates are lower for the BL Lacs. So these can constitute weaker and 
later stages of earlier and stronger FSRQs, provided the cosmic structures evolve so 
as to make less and less mass available for accretion; but this is just what has been 
anticipated above on the basis of the QS evolution alone. 

Second, the link is supported by a closer consideration of the strong evolution shared 
by the FSRQs with the other QSs. For example, Gavaliere & Vittorini (2002, see also refs, 
therein) trace back the accretion episodes that feed the QSs to triggering interactions 
of the host galaxies with their companions in a group. These events destabilize the host 
gas and start inflow toward the nucleus; this develops over a galactic dynamical time 
scales of a few 10“^ Gyr, but some 5 similar cycles are expected per host galaxy since 
z w 2.5. The outcome is exhaustion of the host gas on a time scale of 2-3 Gyr. 

Two consequences are expected. For one, the efficiency of such episodes drops on the 
same scale; so the average QS luminosities are bound to halve on the scale tl « 2 Gyr, to 
be observed as strong LE in the population. Moreover, the frequency of such interactions 
is bound to dwindle as the groups evolve into clusters with lower galaxy densities and 
higher velocity dispersions; this gives rise to positive DE developing over the longer 
scales To ~ several Gyrs set by the hierarchical clustering. These two evolutionary 
components are recognized in the QSO data, see Boyle et al. (2000). 

Thus we expect the powerful FSRQs will fade out over times of some 10“^ Gyr after 
a “last interaction” and a last episode of accretion of mass with associated angular 
momentum. In some 1/2 of the instances this will add constructively to the pre-existing 
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3 , leaving behind a fast spinning hole and a long lived BL Lac. 

Thus the scale for BL Lac births will be close to —To, i-C., the negative of that for 
deaths of bright FSRQs. In the formers’ counts this offsets the slow LE, to the effect 
of flattening them below the Euclidean slope as illustrated in fig. 2. We add that the 
count normalization, i.e., the object number, is bound to be at least a factor 5 below 
the FSRQs’ which undergo repeated activity cycles. 


6. Endpoints 

Note that flattening of the counts is also contributed at faint fluxes by the shape of the 
LF, which is itself flattened at its lower end by the beaming effects (see Urry & Padovani 
1995); but this hardly could swamp a strong evolution of BL Lacs, if it were present. In 
fact, the similarly affected LF of the FSRQs yields counts that do show a steep bright 
section just corresponding to a QS-type intrinsic evolution. 

In the clear absence of the latter, there is scope in focusing on the negative BL Lac 
evolution that is now emerging from the statistical noise in larger and deeper surveys 
with substantial redshift information (Perri et al. 2002). If nailed down at the quantita¬ 
tive levels predicted above, we stress it will provide statistical evidence that transitions 
FSRQs ^ BL Lacs do occur. 

This will clearly mean a key step toward Grand-Unification of the Blazars. It will 
also have implications for the radio sources. If FR I and FR II constitute the parent 
population of the BL Lacs and of the FSRQs, respectively, then similar transitions 
from II to I - once a forbidden proposition, and now one more widely entertained, see 
Padovani & Urry (2001) and refs, therein - will be made acceptable or even likely. 

The other, direct evidence toward Grand-Unification may be provided by finding at 
z ^ 1 the transitional Blazars implied by the genetic connection from FSRQs to BL 
Lacs. We expect that, when caught in the act of switching at regimes m ~ 10“^, such 
objects will show not only intermediate L, but also emission lines and BBB still shining; 
when the latter is subtracted from the continuum in the manner discussed by D’Elia & 
Padovani (2002), the residual non-thermal SED ought to be peaked at iZp beyond the 
optical frequencies. Not many of these objects are to be expected, since by definition they 
will be weaker than the canonical FSRQs while the transition will take times shorter, if 
anything, than the lifetime at the top of to. A number of these objects may have already 
been found, see Sambruna, Chou & Urry (2000); Perlman et al. (2001); Padovani et al. 
( 2002 ). 


Table 1 - The Blazar Main Sequence 



FSRQs 

BL Lacs 

^ CR accelerators 

optical features 
power 
evolution 
top energies 
hole vs. disk 
key parameter 

em. lines, bump 

L < 10^® erg s-i 
strong 

hv ~ 10 GeV 

Lk ^ Ld 
m ^ 1 

no lines & bump 

L < lO"'® erg s“^ 
weak if any 
hv ~ 10 TeV 

Lk Ld 
rh ~ 10“^ 

none 

L < 10*^ erg s-i 
negligible 

Sm ~ 10^“ eV 
vanishing Lk, Ld 
hi 10 “"^ 


Waiting for these two lines of evidence to consolidate, an eye should be kept open on 
a further connection with sources of high energy cosmic rays. This is conceivable and 
even likely if we carry the BMS and the genetic link to their extreme; starting from the 
relation 7m oc one expects values of order 10^° in cases where the total output is 

L ^ lO"*^ erg/s. 

So one may envisage endpoint objects (“dead BL Lacs”) with very low residual to < 
10“^ and possibly in full ADAF conditions, which would feature very faint if any e.m. 
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^ FSRQs 


BL Lacs CR accs.? 


Fig. 3. To illustrate the BMS, and the possible genetic connection from the hyperactive FSRQs 
at high-medium a to the moderate or weak, low a BL Lacs. This may be extended down to 
2 « 0 to including as local relics the objects, e.m. silent but still effective as particle accelerators, 
discussed in §6 


emission. But they would still support strong, nearly unscreened electric fields under 
the widely entertained assumption that magnetic fields of order B ~ 10^ — 10^ G can 
still be held by such vestigial disks. So they can accelerate particles including protons 
up to limiting energies £m ~ ^eBrg {rg/r)°-'^^ ^ 10^° Mg-^® B 3 eV. Such limiting 

energies, long recognized to be accessible to rotating holes, gratifyingly cover the upper 
range of the ultra-high energy cosmic ray spectrum. 

Output levels L ^ 10^^ erg/s as envisaged above can just provide the observed 
UHECR flux; tens of these accelerators could lie within some 50 Mpc and evade in 
the simplest way the GZK cutoff (Boldt & Ghosh 1999). Intergalatic magnetic fields of 
nG strength would blur the geometrical memory of the sources for most except for the 
closest UHECR events. 

We summarize in Table I the main Blazar features that the BMS can explain or pre¬ 
dict. In fig. 3 we schematically illustrate the further steps warranted by negative evolu¬ 
tion and transitional objects: Grand-Unification of Blazars, and the link with UHECR 
accelerators. 
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